Laser-cluster interaction: x-ray production by short laser pulses 
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We investigate the heating of the quasi- free electrons in large rare gas clusters (A'^ exceeding 
10^ atoms) by short laser pulses at moderate intensities (7 ~ lO^^Wcm"^). We identify elastic 
large-angle backscattering of electrons at ionic cores in the presence of a laser field as an efficient 
heating mechanism resembling the Fermi shuttle. Its efficiency as well as the effect of collective 
electron motion, electron-impact ionization and cluster charging, are studied employing a mean- 
field classical transport simulation. Results for the absolute x-ray yields are in surprisingly good 
quantitative agreement with recent experimental results. 
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The interaction of short and ultra-short intense laser 
pulses with clusters has become an important area of 
laser-matter research bridmng the gap between gas-phase 
and solid-state processes The observation of char- 
acteristic x-ray emission from laser-irradiated clusters 
IE ^ IB| suggested its potential as an x-ray source 
through a highly non-linear conversion of IR radiation 
that combines advantages of both solid and gaseous tar- 
gets. Like solids, clusters provide large x-ray yields, yet 
unlike solids they are debris-free, just like gaseous tar- 
gets. Characteristic x-ray emission also provides impor- 
tant time-differential information on the laser-induced 
electronic dynamics on a femtosecond scale. The charge 
state as well as the vacancy distribution of the clus- 
ter ions at the instant of emission can be extracted 
from high-resolution x-ray spectra. Moreover, as the 
vacancy production in deeply bound shells (e.g. K-shell 
in argon or L-shell in xenon) proceeds via impact ion- 
ization by energetic electrons, characteristic x-rays pro- 
vide an "in situ" thermometer of the temperature of the 
heated electron gas in the cluster. Recent experiments 
lEQIl] found an unexpectedly low laser intensity thresh- 
old for the production of x-ray radiation. When irra- 
diating large argon clusters with N > 10^ atoms with 
laser pulses with a short pulse duration of r = 60 fs at 
full width half maximum (FWHM) , the intensity thresh- 
old for the production of characteristic K-x-rays lies at 
Ith = 2.2 • 10^^ Wcm~^. By comparison, the ponderomo- 
tive energy. Up = F^/(4a;^), of a free electron in a laser 
field of this intensity is as low as Up ~ 130 eV, more 
than an order of magnitude below the binding energy 
Ek — 3.2 keV of the K-shell electrons. This observa- 
tion raises puzzling questions as to the efficient heating 
mechanism for electrons in large clusters at such mod- 
erate intensities of very short pulses with ~ 40 optical 
cycles. 

Several theoretical models for intense laser-cluster in- 
teraction have been proposed 0, 0, S IS ^1 0| , none 




FIG. 1: Time evolution of the velocity of an electron in a 
laser field, schematically. A free electron does not experience 
effective velocity gain beyond the quiver velocity (a). If the 
velocity is flipped in synchrony (b), an electron can rapidly 
absorb kinetic energy from the laser field. 



of which appears to provide a satisfactory explanation 
for such rapid acceleration of electrons. A theoretical de- 
scription of intense laser-cluster interaction represents a 
considerable challenge in view of the many-body nature 
of this process. Molecular dynamic simulations 0) S 
are limited to about 1000 atoms, and results obtained for 
small clusters are difficult to scale to larger sizes. The 
recently proposed microscopic particle in cell (MPIC) 
method [lOj reaches clusters of ~ 10** atoms. Larger 
clusters with N > 10^ particles appear still not in reach, 
and quantitative predictions for x-ray emission and inner- 
shell processes have not yet been attempted. We propose 
in this letter an efficient heating mechanism of electrons 
in large clusters that is operational within a few optical 
cycles and at moderate laser intensities. It is based on the 
observation that elastic large-angle scattering of electrons 
at cluster atoms (ions) in the simultaneous presence of a 
laser field provides an efficient route to electron acceler- 
ation. Elastic backscattering at the core potential of the 
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ions can flip the velocity vector of an electron, allowing 
it with non-negligible probability to remain synchronized 
with the alternating laser field vector (Fig. Conse- 
quently, the electron will rapidly accumulate rather than 
lose momentum during subsequent half cycles. Within 
a few optical cycles, electrons can thus be accelerated 
to high kinetic energies well beyond the quiver energy 
Ep — 2Up. This heating mechanism resembles the Fermi 
shuttle acceleration 0, 0| and is also related to the 
lucky-electron model proposed for IR photoemission from 
metallic surfaces The acceleration of a particle by 

successive backscatterings from a moving and a station- 
ary target in the original proposal by Fermi for acceler- 
ation of cosmic particles is here modified such that the 
cluster ions represent stationary targets, while the alter- 
nating force field of the laser plays the role of the mov- 
ing target. A realistic estimate for the efficiency of this 
heating mechanism hinges on a proper description of the 
differential elastic scattering cross sections for electrons, 
dae/d9, into backward angles 9 > 90°, which are deter- 
mined by the non-Coulombic short-ranged potentials of 
the ionic cores. d(Je/d9 was calculated for electron scat- 
tering at argon ions for all relevant charge states between 
q = 1 and q = 16 over a wide range of energies using 
parametrized Hartree-Fock potentials T?, IS 
typically dominated by few low-order partial waves giv- 
ing rise to generalized Ramsauer-Townsend minima [l5j 
and diffraction oscillations [3| (Fig. We assumed 
for simplicity the electronic ground state occupation for 
each q. Extensions to core-excited configurations would 
be straight forward. For the interstitial region a muffin- 
tin potential is employed in order to account for solid- 
state effects 01 . The potential shape in this region has, 
however, no significant infiuence on the cross sections at 
backward angles. The latter exceeds the pure Coulomb 
case by several orders of magnitude for all charge states 
and over a wide range of electron energies (< keV). The 
frequent usage in simulations of unrealistic (softened) 
Coulomb potentials 0, 0, 0| , which grossly underesti- 
mates backscattering, is quite likely one reason why this 
route of electron acceleration has not yet been accounted 
for. Moreover, this process becomes much more impor- 
tant for large clusters, as the mean-free path for elastic 
scattering becomes comparable to the cluster size. The 
important role of realistic core potentials has recently 
also been identified in the quantum analogue of this pro- 
cess, inverse Bremsstrahlung (Up <C hm), for clusters in 
a vacuum ultraviolet (VUV) laser field [l9j . 

A full ab initio simulation for large clusters {N > 10^ 
particles) appears still impractical. In the following we 
present a simplified theoretical description of the elec- 
tronic ensemble that allows to tackle its short-time dy- 
namics (t ~ 60 fs). It employs a generalization of clas- 
sical transport theory (CTT for open systems, in 
which the electronic dynamics is represented by a clas- 
sical phase-space distribution f{r,r,t) whose evolution 
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FIG. 2: (Color online) Left: Differential cross-section distri- 
bution dae{6, E)/dQ, in a.u. for elastic electron scattering at 
ions. The cross section was calculated by partial wave 
analysis for a parameterized Hartree-Fock potential [TtI Il8l|. 
Right: dae{0) / dO. for an electron with fixed kinetic energy 
E — 10 a.u. (solid line). For comparison, the Rutherford 
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is also displayed (dashed line). 



is determined by test-particle discretization, i.e. by solv- 
ing the corresponding Langevin equation for represen- 
tative trajectories. In the present case, the ensemble 
consists of Ne quasi-free electrons, liberated inside the 
cluster after ionization of the cluster atoms, represented 
by Attest = aA^e particles, where the scaling parameter a 
is limited by computational feasibility. Each test particle 
is subject to a Langevin equation (atomic units are used 
unless otherwise stated), 



-FL{t)-VV{r,,t) 



(1) 



with i = 1, . . . Ntest{t)- Eq. n describes a dynamical sys- 
tem open to both particle number variation, Ntest{t), due 
to successive ionization events, and energy exchange with 
the many-particle reservoir (atoms, ions, and electrons) 
as well as with the laser field taken to be linearly polar- 
ized with a temporal envelope 
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Eq. ^ provides a computational starting point for treat- 
ing many-body coUisional correlation effects through 
stochastic forces Fstoa which can be determined either 
from independent ab-initio quantum calculations or ex- 
perimental data [i^l . Forces resulting from static conser- 
vative potentials V inside the cluster , can be included 
as well. Because of the coherent motion and high ion- 
ization density, the present extension of the CTT (Eq. 
^ goes beyond the independent-particle description by 
including dynamic electron-electron and electron-ion in- 
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teractions on a time-dependent mean field level. Ac- 
cordingly, Fraean depends on the entire ensemble of test- 
particle coordinates {ri{t)}. We therefore propagate the 
elements of the ensemble by self-consistently coupling a 
given trajectory to the mean field of other trajectories 
running in parallel by continuously updating the forces. 

Effects of fluctuations on the electronic dynamics can 
be taken into account through stochastic forces which 
are determined from Poissonian random processes. For 
example, electron-ion scattering, electron-impact ioniza- 
tion, and core-hole excitation are determined by the 
probability per unit pathlength for scattering 



dR, 



dx 
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(3) 



controlled by the energy {E) and charge state [q) depen- 
dent integral cross-section for this process, (Tscatt{q, E), 
and the instantaneous ionic target density p{t) of a given 
charge state. Each stochastic scattering process results 
in "jumps" (classical trajectory jumps and jumps in oc- 
cupation) at discrete times. A jump in momentum, Ar, 
signifies elastic scattering determined by the differential 
cross section, dae(q,E,9)/d0, a simultaneous jump in 
test- particle number, ANtest, represents ionizing colli- 
sions, and a simultaneous jump in the number of inner- 
shell vacancies, ANk, results from core-exciting colli- 
sions. The key point is that the necessary input data, 

scatt 

{q,E), can be determined and tabulated indepen- 
dently from the simulation at any desired level of sophis- 
tication. In the present simulation, the electron-impact 
ionization cross sections are determined from a modified 
Lotz formula HHJ cr,(g, E) = A* \n{E/W*)/{EW*) (for 
E > W*), where the empirical parameters A* and W* 
were obtained by a fit to experimental ion-atom collision 
data 0, . The cluster-specific effects of suppression 
of the work function can be incorporated by modifying 
the effective work function W* in the presence of nearby 
ionic cores. 

For the mean field we perform a multipole expan- 
sion keeping only the monopole and dipole terms. The 
monopole term is given by F^°^j,„(r, = {Q{r,t))r /r^ 
where {Q{r,t)) is the instantaneous charge of the sphere 
of radius r resulting from the displacement of the en- 
semble of test particles relative to the ionic background. 
Analogously, the dipole field inside the cluster (r < R{t)) 
is -F^Ln(^) = -P{t) / Ft.{tf , while outside it is that of a 
central dipole. The dipole moment, p is determined by 
p{t) — — ^ Yli where the sum extends over the subset 
of test particles with r; < R{t). 

As the ionic and electronic dynamics proceed on dif- 
ferent time scales, the onset of cluster expansion can be 
taken into account through the parametric variation of 
the radius R{t) of the uniform spherical charge back- 
ground representing the ions of mass M in their time- 
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FIG. 3: Dynamics of an argon cluster with A'' = 2.8 x 10^ 
atoms, irradiated by a laser pulse with wavelength A = 
800 nm, duration r = 60 fs, and peak intensity / = 3.5 x 
10^^ Wcm^^. All quantities are shown as a function of time, 
(a) Mean ionic charge state {q{t)), (b) cluster charge per 
atom Q{t)/N, (c) comparison between plasma frequency LL}p{t) 
(solid line) and laser frequency u (dashed line), (d) effective 
field F^ff inside the cluster including dielectric polarization 
(solid line) and laser field Fl (dashed line). 



dependent mean charge state {q{t)): 

.d^R{t) {q{t)){Q{R,t)) 



M- 



R^t) 



(4) 



We solve (Eq. ^ for a cluster with iV = 2.8 x 10^ ar- 
gon atoms with initial atomic number density p(t = 0) = 
2.66 X 10^^ cm"-^ and initial radius R{0) = 258 a.u. irradi- 
ated by a laser pulse of length (FWHM) r = 60 fs, wave- 
length A = 800nm and peak intensities / > 10^^ Wcm~^. 
The ponderomotive energy of the electrons is of the or- 
der oi Up > 70 eV Ki 2.5 a.u. As the laser field reaches for 
the first time the threshold field for over-barrier ioniza- 
tion Fl{ti) ~ 0.08, the first Ntest{ti) test particles with 
zero velocity randomly distributed over the cluster pro- 
vide initial conditions for the propagation of Eq. In 
the present case Ntest — O-lA^e (or a = 0.1). Contribu- 
tions from tunneling ionization can be included but are 
in the present case negligible. The test-particle number 
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subsequently increases by further ionization events (Fig. 
|3fa)). Impact ionization by electrons meanwhile accel- 
erated is highly effective, making further field ionization 
events unlikely. The mean charge state rapidly increases 
to ^ 7, which most likely still underestimates the ioniza- 
tion efficiency, as non-radiative core-hole relaxation and 
enhanced ionization by suppression of the work function 
by ion proximity are not yet included. The Coulomb ex- 
pansion of the cluster sets in slowly due to the large iner- 
tia of the ions. Even after the laser pulse is switched off 
{t > 2t), the cluster has expanded by less than a factor 
2. Our simulation shows that the charge resulting from 
electrons leaving the cluster (Fig. Efb)) is concentrated 
on the surface of the cluster, the ions in the inside of the 
cluster being well shielded by the quasi- free electrons, in 
agreement with the MPIC simulation After the first 
ionization burst the electronic plasma frequency is given 
by ^2 = Ntestiti)/iaRiti)^) = p{h)ATr/3 ~ 0.016a.u. « 
5oj'^. As electron- impact ionization produces more quasi- 
free electrons, ujp grows rapidly before diminishing again 
as the cluster expansion sets in (Fig. El^c)). During the 
evolution, the effective field inside the cluster consisting 
of both the laser and the polarization field is approxi- 
mately given by 
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where Auj is the Fourier width caused by the temporal 
profile of the pulse (Eq. [2Jl and 7 stands for the damp- 
ing due to scattering events, -fe// is significantly re- 
duced compared to the bare laser field (Fig. Eld)) due 
to the combined effect of collective electron motion and 
electron-impact ionization. A significant resonant en- 
hancement is absent. Ref. 01 suggested laser de- 
phasing heating (LDH), in which the phase shift between 
dipole moment and laser field caused by the macroscopic 
electric field inside the cluster, would lead to a net elec- 
tron energy absorption / JFeffdt, where J stands for 
the electron current. However, evaluating this integral 
using the simulated current and field shows that heating 
by LDH is not sufficient in the present case. 

The efficiency of heating by elastic electron- ion scatter- 
ing is directly refiected in the simulated absolute x-ray 
yields (Fig. 0)). The latter are determined by the num- 
ber of K-shell vacancies created, corrected for the mean 
fiuorescence yield rj taken to be 77 « 0.12 for argon 
with partially filled L-shell but empty M-shell. It should 
be emphasized that the simulation contains no freely ad- 
justable parameter. To compare the simulation results 
to the experiments, an ensemble average over the spatial 
intensity profile of the laser beam, which is Gaussian to a 
good degree of approximation 4], is performed. To quan- 
tify the significance of the Fermi-shuttle acceleration by 
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FIG. 4: Absolute x-ray yield as function of the laser peak 
intensity. Experimental results ;4] (O)) simulation results 
neglecting elastic electron- ion scattering (a) and including 
elastic electron-ion scattering (■). (Lines to guide the eye). 



repeated elastic backscattering, we performed an other- 
wise identical simulation with elastic electron-ion scatter- 
ing switched off. In this case, for / > 1.5 x 10^^ Wcm^^ a 
small fraction of quasi-free electrons gains sufficient en- 
ergy to produce K-shell vacancies. Their mean kinetic 
energy can be estimated from the potential energy of the 
charged- up cluster with charge Q{t) (i.e. the monopole 
term of the mean effective field). However, including elas- 
tic electron-ion scattering drastically increases the x-ray 
yield by a factor 3 to 6. We then find surprisingly close 
agreement with the experimental results (Fig. 01 . 

In summary, we have analyzed the heating of the quasi- 
free electrons in large rare-gas clusters {N 10^ atoms) 
at moderate laser intensities (/ = 10^^ — lO^^Wcm"^). 
We have identified a novel, highly efficient electron heat- 
ing mechanism operative at short times within a few 
optical cycles in terms of elastic large-angle scattering 
resembling the Fermi shuttle. Other processes such as 
heating in a plasma-resonant field are found to be less 
effective. In particular, the polarization of the cluster 
leads to a reduction rather than an enhancement of the 
effective field. While the surprisingly good quantitative 
agreement on an absolute scale with experimental data 
may be, in part, fortuitous, the importance of this route 
to fast electron acceleration appears unambiguously es- 
tablished. 
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